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, Abstract. We investigate linear matter density perturbations in models of structure 

formation with causal seeds. Under the fluid approximation, we obtain the analytic 
solutions using Green-function technique. Some incorrect solutions in the literature 
are corrected here. Based on this, we analytically prove that the matter density per- 
turbations today are independent of the way the causal seeds were compensated into 
J> ■ the background contents of the universe when they were first formed. We also find 

C , that the compensation scale depends not only on the dynamics of the universe, but 

also on the properties of the seeds near the horizon scale. It can be accurately located 
, by employing our Green-function solutions. 
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With the cosmological principle as the basic premise, there are currently two main 
paradigms for structure formation of the universe — inflation [1] and topological 
defects [2]. Although recent observations of the cosmic microwave background 
(CMB) seem to have favored inflation [3], defects can still coexist with it and their 
observational verification will have certain impact to the grand unified theory. In 
the literature the power spectra of models with causal seeds like defects have been 
investigated using the full Einstein-Boltzmann equations. However, the study of the 
phase information of these perturbations still remains difficult due to the limited 
computation power [4]. Although there have been some detailed treatments for 
models with causal seeds [5,6], we shall present a simpler formalism to provide 
not only a physically transparent way to understand the evolution of their density 
perturbations, but also a computationally economical scheme to investigate their 
phase information [7]. This formalism is parallel to those presented in [8] and [9], 
but we give modifications to include the cosmological constant A, as well as some 
other improvements and corrections. 

In a flat Friedmann-Robertson- Walker (FRW) model with an evolving weak 
source field of energy- momentum tensor 0^(x, rj), the full evolution equations of 
linear perturbations can be obtained by considering the stress-energy conservation 
of the fluids and the source, as well as the linearly perturbed Einstein equations [7]. 
With the photon-baryon tight coupling approximation, a closed set of equations for 
the density perturbations in the synchronous gauge are: 
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where R = 3pB/4p r , @+ = ©oo + ®u, a is the scale factor, a dot represents 
the derivative with respect to the conformal time rj, and the subscripts 'c', 'B' 
and V denote cold dark matter (CDM), baryons, and radiation respectively. 1 
By splitting the perturbations into the initial (I) and subsequent (S) parts as 
5/v(x, rj) = 5]y(x, rj) + 5%(~x, rj) where N — c, r, and employing the zero entropy 
fluctuation condition on super-horizon scales as part of the initial condition, we 
solve the above equations for A = in the Fourier space to yield 

^(k, 77) = Gs(k; ri, jft)5 c (k, m ) + G±(k; rj : m )5 c (k, Vi ), (3) 
5%(k, v ) =4ttG f V G^(k; V ,v)^ + (Kv)df,, (4) 

where r]\ is the initial conformal time, and the full expressions of Gf, including 
the baryonic effects, are presented in [7]. We notice that to solve for 5 C we need 
only two Green functions (G% and G%), instead of five as presented in [8], some of 
which are incorrect due to incorrect initial conditions. Figure 1 (left) shows the 
asymptotic behaviors of these two Green functions on the super- horizon (hi) <^ 1) 
and sub-horizon (kfj 3> 1) scales: Tf{k;fj) = lim v / Vcq ^ 00 Gi{k;i],fi)a cq /a, where the 
subscript 'eq' denotes the epoch of radiation-matter density equality. A simple and 
accurate extrapolation scheme can then be used to obtain solutions in the non-flat 
or A 7^ cosmologies [7] . All the solutions are numerically verified to high accuracy. 

One important problem for structure formation with causal seeds is to investigate 
how the source energy was compensated into the radiation and matter background 
when the seeds were formed at 771, and how the resulting perturbations today depend 
on this. First consider the pseudo energy r o = ©00 + (3/87rG) (a/a) 2 (Q C 5 C + 
^r^r) + (a/ a)5 c / AtxG. Since causality requires r o = 3<5 r /4 — 5 C = on super- 
horizon scales, it follows that the initial source energy Goo can be compensated into 
between 5n and 5n with different portions. With our Green-function solutions, it 
can be straightforwardly shown that no matter how @ o was compensated into the 
background at 771, the resulting 5^ and thus 5 C today will be the same. This was 
first numerically observed in [9], and here we can provide an analytic proof. 

Finally we use our Green-function solutions to study the compensation mecha- 
nism and the scale on which it operates. First it can be shown that [7] r o(k, 770) = 



(i-r(fc))eoo(k,ijb)+j; 



T'(k- r))(a(rj)/a(rj))e + (k, rj) + T(A;)0 O o(k, fj) 



drj, where 



T(k) is the standard CDM transfer function in the inflationary models, and 



) Here we have ignored neutrinos, whose effects on the current study is negligible. 
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FIGURE 1. The asymptotic behaviors of the CDM Green- function solutions on the su- 
per-horizon and sub- horizon scales (left), as well as on all scales (right). 



T'(k;fj) = Ql(k;7] ,fi)/Ql Q (k;r] ,fi) (see Figure 1 right). On super-horizon scales, 
it is clear that T(k) is unity by definition so that only the integral survives, and 
that the quantity inside the square brackets is nothing but © i,j(k, fj) due to the 
source stress-energy conservation and the fact T = X". Since Qq^ has a k 4 fall-off 
power spectrum due to causality, it follows immediately that the pseudo-energy t oq 
also has a A; 4 -decay power spectrum outside the horizon. On sub-horizon scales, on 
the other hand, although (1 — T(k)) is approximately unity, the usual sub-horizon 
power-law decay in O o (as in the case of topological defects) will make the first 
term negligible, while the quantity inside the integral is no longer simply ©oi,j- As 
a result, we see that the compensation scale, above which the power of both S c and 
Too decays determined not only by the functions T and T", but also by the 

properties of the source near the horizon scale. Once the detailed behavior of the 
source near the horizon scale is known, we can accurately locate the compensation 
scale using our formalism [7]. We acknowledge the support from NSF KDI Grant 
(9872979) and NASA LTSA Grant (NAG5-6552). 
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